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Present knowledge of the function of materials is largely based on studies (experimental and 
theoretical) that are performed at low temperatures and ultra-low pressures. However, the majority 
of everyday applications, like e.g. catalysis, operate at atmospheric pressures and temperatures at or 
higher than 300 K. Here we employ ab initio, atomistic thermodynamics to construct a phase diagram 
of surface structures in the (T, p)-space from ultra- high vacuum to technically-relevant pressures 
and temperatures. We emphasize the value of such phase diagrams as well as the importance of the 
reaction kinetics that may be crucial e.g. close to phase boundaries. 
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The main prerequisites for reaching a microscopic un- 
derstanding of heterogeneous catalysis are the identifi- 
cation of the composition and geometry of the catalyst' 
surface and the determination of the various chemical re- 
actions that take place under realistic conditions. Unfor- 
tunately, most Surface Science experimental techniques 
are difficult if not impossible to use at the pressures (of 
the order of one atmosphere) and temperatures (often 
higher than 300 K) that are typically applied in steady- 
state catalysis. Therefore, what is considered to be im- 
portant elementary processes at the catalyst' surface (e.g. 
dissociation, diffusion, and chemical reactions) has usu- 
ally been concluded from chemical intuition and exten- 
sive knowledge from ultra-high vacuum (UHV) experi- 
ments. However, several studies revealed the danger of 
this approach (see e.g. Ref. and references therein), 
and the need to reliably bridge the temperature and pres- 
sure gap between UHV and "real life" is probably the 
main challenge in modern Surface Science. 

In ab initio theory the consideration of high temper- 
ature and high pressure can be achieved by explicitly 
taking into account the surrounding gas phase in terms 
of ll ab initio, atomistic thermodynamics" (cf. e.g. Refs. 

||, f|, ||, [|). This is also an appropriate (first) ap- 
proach to steady-state catalysis, which is often run close 
to thermodynamic equilibrium (or a constrained equilib- 
rium) to prevent catalyst degradation. We will qualify 
this statement below when analyzing our results. In the 
following we show how the combination of thermodynam- 
ics and density-functional theory (DFT) can be applied 
to obtain the lowest-energy surface structures in a (con- 
strained) equilibrium with the surrounding gas phase, 
thus enabling us to construct a (T, p)-diagram of the sta- 
bility regions (or metastability regions) of different sur- 
face phases. We use the example of CO oxidation over a 
RuO2(110) model catalyst to illustrate the concepts and 
conclusions that can be derived from such an ab initio 
surface phase diagram. 

At given temperature, T, and partial pressures, {pi}, 
the stable surface structure is the one that has the 



lowest surface energy, j(T,{pi}) = l/A[G(T,{pi}) - 
Y^iNiHiiT^i)}. Here, G(T,{ Pl }) is the Gibbs free en- 
ergy of the finite crystal, JV» and Hi(T,pi) are number 
and chemical potential of the species of the ith type, and 
7(T, {pi}) is measured in energy per unit area by dividing 
through the surface area, A. The surface free energy is 
thus a function of the chemical potentials of all species in 
the system, e.g. in the present application to CO oxida- 
tion over RuO2(I10) we have three fj,i(T,pi) correspond- 
ing to the three chemical elements: Ru, C, and O. Con- 
sidering that the surface is not only in equilibrium with 
the gas phase, but also in equilibrium with the underly- 
ing metal oxide, implies that the chemical potentials of 
Ru and O are not independent variables but tied by the 
Gibbs free energy of the bulk oxide. As a consequence, 
the surface free energy depends only on two chemical 
potentials, which we choose to be those of O and CO. 
Because chemical potentials are directly related to tem- 
perature and partial pressure, a comparison to the whole 
range of experimentally accessible gas phase conditions is 
achieved. [|J Here, we note that the above description 
refers to a constrained thermodynamic equilibrium only 
with the reactants, because it assumes that the O2 and 
CO gas phase are independent reservoirs, and the direct 
(non catalytic) formation of CO2 in the gas phase can be 
ignored because of its negligible rate. 

As shown previously fi, for Ru02 it is a good ap- 
proximation to assume that the vibrational contribu- 
tions to the Gibbs free energy, G(T, {pi}), nearly equal 
the corresponding terms of the bulk oxide. We will 
therefore replace these quantities by the correspond- 
ing total energies, which are calculated using DFT. 
We use the generalized gradient approximation (GGA) 
for the exchange-correlation functional ||, and the all- 
electron Kohn-Sham equations are solved by employ- 
ing the full-potential linear augmented plane wave (FP- 
LAPW) method p0| p|. The high accuracy of these cal- 
culations (cf. Ref. M) implies that the numerical accu- 
racy of the relative j(T, {pi}) is better than ±5 meV/A 2 . 
We also note that using the local-density approxima- 
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tion (LDA) as exchange-correlation functional gives es- 
sentially the same result for the surface phase diagram 
(cf. Fig. ||). Thus, the DFT accuracy for this system's 
surface phase diagram is very high, partly due to the fact 
that it follows only from energy differences. 

Ruthenium is well known to exhibit an exception- 
ally high activity for the catalytic oxidation of CO, 
with turnover rates exceeding those of the more com- 
monly employed materials like Rh, Pd, or Pt.JlJ, [ll| 
For Ru(0001) this high activity could be related to 
RuO2(110) oxide patches that are formed in the reac- 
tive environment. [114 , 15, lq] The atomic arrangement at 



RuO2(110) can be described by starting from the stoi- 
chiometric surface (see the orange-framed inset in Fig. 
0, labeled O br /-). There are two different oxygen atoms 
at this surface: threefold-bonded O (as is an O atom in 
bulk RUO2) and twofold bridge-bonded O br . Likewise, 
there are two different surface Ru atoms, one of which 
is bonded to six oxygen neighbors (as is a Ru atom in 
bulk R11O2), and one is bonded to only five O neigh- 
bors. The latter represents a "coordinatively unsatu- 
rated site" and is called Ru cus . The bulk stacking se- 
quence would be continued by adsorbing additional oxy- 
gen atoms on top of the Ru cus atoms. This gives the 
red-framed oxygen-rich geometry, O br /O cus . The adsorp- 
tion energy for a single O atom on the Ru cus atom of 
the O br /— surface is calculated to be exothermic by 1.2 
eV/adatom (at /i = eV)@- And thc 

energy of an 

oxygen vacancy in the O br /— surface is calculated to be 
endothermic by 2.5 eV/vacancy (removing a O br atom) 
and « 3.5 eV/vacancy (removing a 3f atom) at /to = 
eV[fL7|. Removing all bridge-bonded oxygen atoms from 
the O br /— surface gives the black-framed, oxygen-poor 
geometry — /— in Fig. [IJ which represents the third pos- 
sible (1 x 1) termination of RuO2(110). 

The /io scale can be converted into a pressure scale 
at any chosen temperature,^, [?[] and this is noted at the 
top of Fig. [l] for T = 600 K, a typical annealing temper- 
ature used experimentally in this system Jl5|, [ll| . From 
the computed surface free energies of the aforementioned 
three RuO2(110) terminations, we obtain that even at 
the lowest possible O2 pressures, O atoms will at least 
occupy the bridge sites, leading to the stoichiometric sur- 
face termination (the orange line tagged as O br /— in Fig. 
|l|). At even lower ^o(T,po 2 ) than the O-poor boundary 
in Fig. [l], Ru02 will decompose into Ru metal. || 

We now consider the presence of CO as a second 
species, discussing the surface structure of Ru02 in the 
constrained equilibrium with a CO and O2 gas phase. For 
this we computed what we believe are all possible (lxl) 
surface phases including O and CO at the surface. |l8| 
Figure [l] shows the corresponding surface energies in thc 
CO-rich limit for /ico of three CO containing structures 
that will play a role in the later discussion. In this 
CO-rich limit and if the O chemical potential is below 
— 1.2 eV, CO is transformed into graphite and 1/2 O2. 
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FIG. 1: Surface free energies, y(T, {pi}), in the experimen- 
tally accessible range of the oxygen chemical potential.^, [jj] 
Shown are three RuO2(110) terminations (dashed lines and 
black, red, and orange insets), and the three most stable CO 
containing surface phases for the CO-rich limit (solid lines 
and gray, green, and blue insets). The labels note whether 
bridge or cus sites are occupied by O or CO, or are empty 
(— ). Ru atoms are drawn as light, large, O atoms as dark, 
small, and C atoms as white, small spheres. 



This is why the CO related (gray, blue, and green) lines 
in Fig. [I] exhibit a kink at /to = —1-2 eV. As appar- 
ent from Fig. |l| a phase with CO occupying all bridge 
sites (the gray line tagged as CO br /— ) is barely more sta- 
ble than the hitherto discussed pure O-terminations; and 
occupation of also the cus sites (CO br /CO cus ) leads to a 
(meta)stable phase only at low O chemical potential. To- 
wards higher /^o we find a third relevant geometry, that 
exhibits O atoms at the bridge sites and CO at the cus 
sites, (O br /CO cus ). 

If we now in addition to /to investigate the depen- 
dence on /ico j Fig- turns into a complicated three- 
dimensional plot. We therefore present our results by 
showing only the lowest-energy surface structures (con- 
cluded from many calculations of the type shown in Fig. 
[I] for any value of /ico)- Figure [| displays this central 
result of our study from which the stability regions of 
all surface phases in (constrained) equilibrium with an 
environment formed by O2 and CO can be seen. At 
very low CO chemical potential corresponding to very 
low CO pressure in the gas phase, we recover our results 
obtained for RuO2(110) just in equilibrium with an O2 
environment (dashed lines in Fig. [I]): At low O chemical 
potential we find the stoichiometric termination with O 
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FIG. 2: Regions of the lowest-energy structures in (po, /-tco)- 
space (see Fig. [l] for the nomenclature). The additional axes 
give the corresponding pressure scales at T = 300 K and 
600 K. In the blue-hatched region gas phase CO is trans- 
formed into graphite. Regions that are particularly strongly 
affected by kinetics are marked by white hatching (see text). 



at bridge sites, while towards higher O2 pressures the O- 
rich termination with O additionally located at the cus 
sites becomes more stable. Increasing the CO content in 
the gas phase, CO is first bound at the cus sites. This 
is easier at low O chemical potentials, where these sites 
are free. At higher /io oxygen atoms also compete for 
cus sites, which is why the O br /CO cus phase becomes 
only more stable than the O br /O cus phase at progres- 
sively higher CO chemical potentials, cf. Fig. [| Finally, 
at a high CO and low O chemical potential in the upper- 
left corner of Fig. ||, we find that CO is also stabilized at 
the bridge sites (instead of O), leading to a completely 
CO covered surface (CO br /CO cus ). 

The knowledge of such a phase diagram (Fig. ^) is 
a first, important step towards an understanding where 
in (T, {pi})-space catalysis may be most efficient, but 
in parallel it also allows to systematically discuss under 
which gas phase conditions our constrained thermody- 
namic equilibrium approach may break down: We re- 
call that the imposed constraints are that (i) in the gas 
phase the two molecules will not react to CO2, which 
is appropriate because the free-energy barrier is signifi- 
cant. Thus, CO2 is only formed at the catalyst surface 
from where it will desorb immediately. Furthermore, it 
is assumed that (n) the CO2 formation at the surface is 
slower than other adsorption and desorption processes, 
so that the surface can maintain its equilibrium with the 
reactant gas phase. Finally, (iii) the gas phase should 



not affect the stability of the RA1O2 bulk structure. 

With respect to the latter point, we note that for 
/^co > MO + A, Ru02 will decompose into Ru metal 
(A th = 0.06 eV, see the dotted white line in Fig. 0; and 
A ox p = 0.26 eV@). Thus, a noticeable part of Fig. | 
(the blue and most of the upper-left green region) refers 
to metastable situations and will not prevail for long un- 
der realistic conditions. As for the remaining part of the 
phase diagram, it is clear from the pressure scales in Fig. 
H that the stoichiometric surface (orange region, denoted 
as O br /— ), that had been observed in UHV experiments, 
is not very important under high-pressure, catalytically 
relevant conditions. Obviously, for oxidation catalysis it 
is important that both reactants (here CO and O) are 
adsorbed at the surface. In this respect, the green region 
(O br /CO cus ) appears to be the most relevant one at first 
glance, and we have to ask if our remaining assumption 
(ii) of the constrained thermodynamic approach is ap- 
propriate here: If the reaction CO cus + O br — > CO2 were 
faster than the dissociative adsorption of O2, adsorbed 
CO would rapidly eat away bridge-bonded oxygen atoms, 
and surface kinetic effects would render the oxygen con- 
centration much lower than the full monolayer that is 
behind the atomic structure of the green region. Fortu- 
nately, for this region the CO2 formation energy barrier 
is calculated to be noticeable (1.2 eV, in agreement with 
Ref. [|l9|), and the energy to create vacancies in the O- 
bridge layer is high (see above). Correspondingly, the 
thermodynamic approach seems to be valid here. 

The question remains where in (T, p)-space do we ex- 
pect the Ru02 catalysis to be most efficient? In this 
context we note that a so-called stable phase is not sta- 
ble on an atomistic scale, but represents an average over 
many processes such as dissociation, adsorption, diffu- 
sion, association and desorption. As all these processes 
and their interplay |lj are of crucial importance for catal- 
ysis, regions in (T, p)-space where such fluctuations are 
particularly pronounced can be expected to be most im- 
portant. This is the case close to the boundaries between 
different phases: At finite temperatures, the transition 
into a neighboring phase occurs not abruptly in (po, 
A*Co)-space, but over a pressure range in which the other 
phase gradually becomes more populated. The resulting 
phase coexistence at the catalyst surface may then lead 
to a significantly enhanced dynamics, in which even ad- 
ditional reaction mechanisms can take place and/or (dy- 
namic) domain pattern formation may occur. Assuming 
a canonic distribution we estimated in Fig. || the region 
on both sides of the boundaries in which the respective 
other phase is present at least at a 10% concentration at 
room temperature (white areas). From the given pres- 
sure scales we see that of particular interest for catalysis 
is the green/red boundary between the O br /CO cus and 
O br /O cus phases, where oxygen atoms and CO molecules 
compete for the same site, i.e. the cus site. 

In fact, in this phase-transition region and in the neigh- 
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boring hatched part of the red phase in Fig. assump- 
tion (it) of our constrained equilibrium theory is not valid 
and the oxygen coverage will presumably be noticeably 
lower than what the present approach suggests: Under 
these conditions adsorbed CO will react with O cus , and 
the filling of empty sites with CO will be fast, whereas the 
filling with oxygen atoms will be slow. In this region high 
catalytic activity is expected and we note that here the 
coverage and structure (i.e., the very dynamic behavior) 
must be modeled by statistical mechanics. However, this 
is not the point of the present paper that aims to iden- 
tify the catalytically active region in (T, {pj})-space, and 
the neighboring stable phases. The exceptionally high 
turnover rates over working Ru catalysts were measured 
for conditions, where both gas phase species are present 
at ambient pressures and about equal partial pressures 
fl2| , [l3[ pp| , i.e. for conditions which lie very close to the 
red-green boundary in Fig. ^, in accordance with our 
conclusions. 

Finally, Fig. || also shows that the pressure gap can 
be bridged when sufficient information about the stabil- 
ity regions of the various surface phases is available, i.e. 
when it is assured that one stays within one phase region 
or along one particular phase boundary. A decrease of 
ambient gas phase pressures over several orders of magni- 
tude maintaining an about equal partial pressure ratio of 
O2 and CO will e.g. still result in conditions rather close 
to the green-red boundary in Fig. |^, and similar reaction 
rates have indeed recently been noticed by Wang et al. 
pl[ when comparing corresponding UHV steady-state ki- 
netic data with high pressure experiments of Zang and 
Kisch p^ j. Yet, from Fig. ||it is also clear, that without 
knowledge of a (T, {pi}) -phase diagram a naive bridging 
of the pressure gap by simply maintaining an arbitrary 
constant partial pressure ratio of the reactants can easily 
lead to crossings to other phase regions and in turn to 
incomparable results. 

In conclusion we computed the phase diagram of sur- 
face structures of Ru02(f 10) in constrained equilibrium 
with a gas phase of O2 and CO. Depending on the tem- 
perature and partial pressures of these reactants, sev- 
eral stable and/or metastable surface phases are found. 
The presented calculations are of highest quality (at 
this time), but the employed GGA is an approximation, 
and the resulting uncertainty may well be of the order 
±100 K. Fortunately this should affect our conclusions 
only little. In particular, our study shows that if one 
stays within the region of one phase or along one partic- 
ular phase boundary, a reliable bridging of the pressure 
gap is possible. While we emphasize the importance of 
such phase diagrams (Fig. ||), we also note that under- 
standing of the full catalytic cycle requires a kinetic mod- 
eling, but for a system as complex as the present one, this 
is not yet possible. 

Regions in (T, {pi})-space that exhibit enhanced ther- 



mal fluctuations, i.e. where the dynamics of atomistic 
processes is particularly strong, are particularly interest- 
ing for the function of surfaces under realistic conditions 
- not only for the example discussed in this paper. The 
described approach identifies such regions close to bound- 
aries in the computed surface phase diagram, where we 
correspondingly expect a particularly high catalytic ac- 
tivity. For our example of the CO oxidation reaction 
over RuO2(110) this conclusion is indeed consistent with 
existing data |1, ll| M. 
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